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ABSTRACT 


ANECHOIC  CHAMBER  SIMULATION  AND  VERIFICATION 

by 

ROBERT  M.  TAYLOR 

Master  of  Science  in  Electrical  Engineering 


This  project  considers  an  electromagnetic  fieid  simulation  of  an  anechoic  chamber,  electromagnetic 
field  simulation  with  verification  measurements.  The  simulation  is  a  Geometric  Optics  (Ray  Tracing) 
mathematical  model  of  the  direct  path  between  two  antennas  and  interfering  scattering.  There  are 
two  separate  models  due  to  the  frequency  dependent  nature  of  the  pyramidal  radar  absorbing 
material  (RAM).  The  model  for  the  frequency  range  of  30  to  500  MHz  was  used  to  characterize  the 
lossy  specular  scattering.  The  specular  scattering  was  modeled  as  a  lossy,  tapered,  TEM 
transmission  line  in  an  inhomogeneous  anisotropic  tensor  material.  The  frequency  range  from  500 
MHz  to  18  GHz  was  characterized  by  dominant  tip  diffraction  of  RAM  patches  and  the  model  made 
use  of  a  Uniform  Theory  of  Diffraction  code  for  a  dielectric,  corner.  The  measurements  and 
simulations  were  based  on  an  azimuthal  cylindrical  scan.  Diagnostic  measurements  were  also 
performed  by  a  cylindrical  scan  of  a  directional  horn  antenna  to  locate  scattering  sources  in  the 
chamber.  A  cylindrical  wave,  modal  expansion  of  the  diagnostic  data  which  included  a  one 
dimensional  Fast  Fourier  Transform  with  Hankel  function  expansions. 
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ANECHOIC  CHAMBER  SIMULATION  AND  VERIFICATION 


INTRODUCTION 

The  performance  of  an  anechoic  chamber  can  be  limited  by  the  imperfections  in  the  radar 
absorbing  materials  (RAM)  or  other  chamber  scattering  sources.  Measurements  to  characterize  the 
imperfections  are  conventionally  in  quantities  of  amplitude  or  phase  ripple,  amplitude  taper  or  site 
attenuation  deviation  (SAD)  within  a  test  zone  or  quiet  zone.  Diagnostic  measurements  such  as 
free-space  standing  wave  ratio  (SWR)  measurements,  site  attenuation  comparisons,  time-domain 
reflectometry,  pattern  comparisons,  plane  wave  spectral  (PWS),  and  inverse  synthetic  aperture  radar 
(ISAR)  imaging,  can  be  used  to  locate  interfering  sources  and  corrections  can  be  made. 

It  is  useful  to  predict  the  electromagnetic  field  characteristics  of  the  chamber  with  mathematical 
models  to  aid  in  the  design  of  the  chamber  or  the  strategic  placement  of  RAM.  The  Geometric 
Optics  (GO)  or  ray  tracing  technique  has  been  applied  to  tapered  chambers  [1],  rectangular 
chambers  [2]-[8].  low  frequency,  site  attenuation  chambers  [9],  and  to  indoor  compact  ranges  [10]. 
A  modification  to  the  conventional  ray  tracing,  infinitely  thin,  ray  model  to  a  reflection  cone  model 
has  been  suggested  [il]  to  include  curved  surfaces  in  a  chamber.  Ray  tracing  has  been  combined 
with  Uniform  Theory  of  Diffraction  (UTD)  for  a  dielectric  corner  in  a  compact  range  [12],  [13].  A 
finite  element  formulation  has  been  suggested  for  anechoic  chamber  analysis  [14]  though  the  initial 
study  was  in  waveguide.  The  method  of  moments  has  been  applied  to  the  analysis  of  four-wire 
anechoic  chambers  [15]  and  [16].  The  finite-difference  time  domain  technique  has  been  used  to 
study  four-wire  anechoic  chambers  [16]  and  rectangular  chambers  [3]  and  [17],  and  a  finite- 
difference  frequency  domain  technique  has  been  useful  in  analyzing  low  frequency  shielded  semi- 
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anechoic  chambers  [18].  A  Transverse  Electromagnetic  (TEM)  wave  model  has  also  been  found 
useful  in  analyzing  four-wire  chambers  [19].  Time  domain  predictions  have  been  used  to  improve 
a  compact  range  performance  [20].  Finite-range  gain  equations  have  been  used  to  evaluate 
chamber  range  antenna  effectiveness  [21].  A  model  management  system  has  been  developed  [22] 
for  Electronic  Warfare  (EW)  simulation  type  chambers  to  include  not  only  the  multipath  interference 
of  the  chamber  walls,  but  also  includes  effects  as  System  Under  Test  (SUT)  interaction  with  the 
chamber,  and  pulse  transmission  characteristics. 

A  key  component  of  the  chamber  simulation  was  the  RAM  model  for  which  various  methods  have 
been  used  to  predict  the  RAM  characteristics.  Some  models  specifically  generated  for  pyramidal 
RAM  are  tapered  transmission  line  models  with  increasing  attenuation  toward  a  conducting  wall 
[23]-[25],  UTD  for  corner  diffraction  [12]  and  [13],  propagation  in  uniform  and  lossy  dielectric  layers. 
Boundary  Element  Method  (BEM),  radiating  structures,  periodic  method  of  moments  (PMM),  finite 
element,  and  statistical  methods,  of  which  the  first  two  techniques  were  used  extensively  for  this 
project.  The  following  paragraph  gives  a  brief  overview  of  the  highlights  of  the  project. 

This  project  was  a  study  of  an  anechoic  chamber,  computer  simulation  for  predicting  the 
electromagnetic  field  characteristics  within  a  test  or  quiet  zone  (figure  1).  Samples  of  the  field  were 
collected  by  a  cylindrical  scan  and  plotted  in  amplitude.  The  effects  of  amplitude  taper,  amplitude 
ripple  and  phase  were  included.  Ray  tracing  was  used  to  simulate  a  propagating  wave  from  a 
stationary  range  antenna  to  a  receive  probe  antenna  over  a  direct  path,  the  specular  paths  from  the 
walls,  floor  or  ceiling  forward  of  the  range  antenna,  and  paths  to  pyramidal  RAM  tips.  The 
simulation  was  divided  into  low  and  high  frequency  models  because  of  the  frequency  dependant 
nature  of  the  RAM  characteristics.  The  low  frequency  model  was  most  effective  in  the  frequency 
range  of  30  to  500  MHz  and  the  high  frequency  model  was  used  from  0.5  to  18.0  GHz.  The  low 
frequency  model  is  similar  to  that  used  by  Kuester,  et.  al.  [9],  in  that  a  ray  tracing  model  was  used 
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and  a  similar  RAM  model  was  used.  However,  a  uniform  layer  at  the  base  of  the  RAM  and  a 
concrete  dielectric  layer  under  the  floor  RAM  was  added  in  this  project.  Other  differences  include 
the  addition  of  the  twisted  pyramid  RAM  on  the  floor,  and  the  output  data  format  of  a  spatial, 
cylindrical  scan  instead  of  the  site  attenuation  frequency  scan  in  [9].  The  high  frequency  model  was 
similar  to  the  chamber  model  in  [12]  and  [13]  in  that  ray  tracing  was  used  and  the  bistatic  reflection 
derivation  is  the  same.  The  bistatic  pyramidal  RAM  tip  diffraction,  reflection  coefficient  subroutine 
in  the  file  TIPBIST.FOR’  from  CHAMBER’  was  linked  and  adapted  to  the  high  frequency  model. 
The  high  frequency  chamber  simulation  differs  from  CHAMBER  as  follows:  The  was  developed  for 
a  small  compact  range  model  whereas  this  project  was  adapted  to  a  large  rectangular  chamber. 
Time  gating  was  not  used  in  this  project  as  in  CHAMBER.  The  chamber  was  not  symmetric  in  all 
cases  in  this  project  whereas  symmetry  was  used  to  reduce  the  number  of  computations  in 
CHAMBER.  The  application  configuration  was  bistatic  in  this  project  instead  of  monostatic.  The 
RAM  patch  configuration  differs  slightly,  also. 

A  verification  of  the  effectiveness  of  the  chamber  simulation  was  conducted  by  comparing  the 
simulation  with  measurements  performed  in  the  chamber  by  a  cylindrical  field  probe  scan.  The 
verification  was  conducted  at  300  MHz  for  the  low  frequency  model  and  10  GHz  for  the  high 
frequency  model.  Differences  in  the  simulation  and  measurements  were  resolved  by  processing  the 
measured  data  by  a  cylindrical  mode  expansion  to  determine  the  angle  of  chamber  scattering.  The 
cylindrical  mode  expansion  method  was  similar  to  that  used  by  Bennet,  et.al.  [26]  which  was  applied 
to  data  collected  on  an  outdoor  range. 

The  following  section  will  describe  the  anechoic  chamber  geometry.  The  differences  between  the 
low  and  high  frequency  model  will  first  be  explained,  followed  by  the  RAM  geometry,  and  then  the 


’Joseph,  P.,  "A  UTD  Scattering  of  Pyramidal  Absorber  for  Design  of  Compact 
Range  Chambers",  Ph.D.  Dissertation,  OSU,  1988. 


3 


wall,  ceiling  and  floor  geometry. 


1.  CHAMBER  MODEL  GEOMETRY 


The  computer  simulation  project  was  developed  to  support  a  rectangular  anechoic  chamber  with 
a  80  ft.  diameter  turntable  located  on  the  floor  with  the  main  coordinate  system  positioned  at  the 
center  of  the  chamber  turntable.  The  application  of  the  simulation  was  for  a  large  rectangular 
chamber  (figure  2)  with  internal  dimensions  of  70  feet  high,  250  feet  wide  and  265  feet  long  with  an 
80  foot  diameter  turntable.  Pyramidal  RAM  covers  all  of  the  wall,  floor  and  ceiling  of  the  chamber 
with  18  inch  long  pyramids  used  in  all  cases  except  for  the  floor  which  was  48  inch  long.  The  east 
and  west  walls  of  the  chamber  was  located  on  the  y  =  -125  and  y  =  125  ft.  planes,  respectively. 
The  floor  and  ceiling  are  located  on  the  2  =  0  and  z  =  70  ft.  planes,  respectively.  The  back  wall 
was  located  on  the  x  =  -109.2  ft.  plane. 

The  probe  antenna  was  located  on  the  turntable  (figure  3)  and  was  scanned  in  the  +0t  direction. 
Figure  4  defines  R*,  R<,.  and  r.  The  range  antenna  location  was  given  by  equation  (1).  A  cylindrical 
scan  was  formed  when  the  probe  location  radial  dimension  was  held  constant  as  the  height  was 
varied.  The  direct  path  between  the  two  antennas  is  given  by  equation  (2).  The  cylindrical  scan  is 
given  in  equation  (3).' 


X  *  Ay  9  ^  Z 

(1) 

II 

1 

(2) 

<})j.  +  lr|sin  4)j,  ^  +  h  z 

(3) 
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Figure  1  Anechoic  Chamber  Configuration  Showing  Multipath 
Effects  in  the  Quiet  Zone 


Figure  2  Wall  Geometry 
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Source 

Antenna 


Azinuth  Rotation 


Figure  3  Cylindrical  Scan 


Figure  4  Cylindrical  Scan  Geometry 
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A.  I  nw  Ffeouenrv  and  High  Frequency  Model  Geometry  Contrast  The  wall,  floor  and  ceiling 
multipath  interaction  configuration  will  differ  in  the  modeling  of  the  low  frequency  characteristics  and 
the  high  frequency  characteristics.  The  electromagnetic  differences  will  be  the  subject  of  a  later 
section  of  this  report.  The  vector  geometry  differences  between  the  low  and  high  frequency  model 
will  be  discussed  here  and  then  the  application  to  each  wall,  floor  or  ceiling  will  be  described. 


The  low  frequency  model  was  an  approximation  of  the  attenuation  of  the  specular  reflection  (figure 
5)  since  near  the  low  end  of  the  frequency  spectrum,  30  MHz,  the  length  of  the  RAM  was  much  less 
than  a  wavelength.  This  implies  there  was  one  vector  to  each  wall,  floor  or  ceiling  per  scan  point 
from  the  range  antenna  to  the  coordinates  of  the  wall,  floor  or  ceiling,  and  one  vector  from  the 
conducting  wall,  floor  or  ceiling  to  the  scanning  probe.  The  method  of  images  [27,28]  (figure  6)  will 
facilitate  the  configuration  so  that  the  two  vectors  are  combined  into  one  vector  from  the  range 
antenna  image  to  the  probe  antenna.  Since  the  RAM  was  modeled  as  a  bistatic  reflection 
coefficient,  by  Snell’s  law  of  reflections  [27]  (0^  =  0J,  the  reflection  and  material  refraction  are 
symmetric,  and  the  image  ray  vectors  will  trace  a  straight  line  from  image  to  probe.  The  E-Field 
formulation  of  the  multipath  will  include  a  consideration  for  the  effects  of  the  bistatic  reflection 
coefficient,  image  boundary  conditions,  and  spherical  wave  spreading.  The  incident  angle  between 
the  negative  of  the  incident  wave  vector  and  the  normal  vector  from  the  coordinates  of  the  walls, 
floor  or  ceiling  was  computed  by  the  general  equation  in  (4), 

•  A'B  =  |a||.B|cos  0  =  +  A^y  +  A^g 

where  A  and  B  are  arbitrary  vectors  and  O  is  the  angle  between  the  vectors.  The  incident  angle 
was  then  introduced  into  the  low  frequency  model. 

The  high  frequency  model  was  concentrated  entirely  to  the  tips  of  the  pyramidal  RAM  (figure  7) 
where  a  non-specular  scattering  is  dominant  to  the  specular  scattering.  The  specular  component 
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of  the  field  was  shown  by  Joseph  in  [13]  to  be  attenuated  enough  to  be  neglected  at  high 
frequencies. 

It  was  also  shown  by  Joseph  in  [13]  that  each  tip  of  the  RAM  will  contribute  to  the  total  scattering 
in  a  rarxJom  fashion  so  that  the  phase  information  of  each  tip  scattering  will  be  lost.  He  defined  a 
patch  of  pyramids  and  the  scattering  component  of  each  tip  was  summed  for  a  total  scattering 
component  for  each  wall,  floor  and  ceiling  of  the  chamber.  The  concept  of  the  random  nature  of 
the  phase  was  used  in  this  project  when  sufficient  computer  memory  could  not  be  allocated  to 
include  the  entire  wall  of  the  chamber.  The  computer  simulation  of  the  scan  was  designed  to 
accommodate  large  scan  radii,  which  required  large  areas  of  the  chamber  to  be  modeled. 
Therefore,  in  lieu  of  defining  an  exact  location  of  each  pyramid  tip,  a  new  patch  was  defined  for 
each  scan  angle.  The  new  patch  was  centered  at  the  specular  location  of  the  scan.  Each  patch 
of  RAM  surrounding  the  specular  point  was  made  as  large  as  the  mainframe  computer  would  allow. 
At  each  scan  angle,  the  E-field  phasor  from  each  tip  was  then  added.  From  the  random  nature  of 
the  phase  from  each  tip,  it  follows  that  the  error  introduced  by  this  method  was  kept  to  a  minimum. 

The  procedure  for  defining  a  typical  patch  (figure  8)  is  as  follows:  A  vector  from  the  range  antenna 
was  considered  to  be  incident  at  a  specular  location  at  the  tip  of  a  pyramid  which  was  at  the  center 
of  the  patch.  The  size  of  the  M  x  N  patch  was  specified  in  odd  numbers.  The  surrounding 
pyramids  in  the  patch  was  oriented  around  the  center  patch  by  the  relationship  in  equations  (5)  and 
(6). 

‘<3  a-(j-l)  (5) 

Yj,  =  -  -a  +  a-(k-l)  (6) 
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Figure  8  RAM  Patch  Top  View 
Showing  Local  Coordinate  System  of 
Tip  5 


Figure  7  Dominant  Tip  Diffraction  for 
High  Frequency  Model 
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where  i  =  1  to  M,  j  =  1  to  N,  and  a  is  the  spacing  between  the  RAM  tips  for  the  patch.  A  cartesian 
coordinate  system  was  defined  for  each  RAM  tip  with  the  z-axis  along  a  unit  normal  vector  from  the 
coordinates  of  the  wall,  floor  or  ceiling.  The  x  and  y  axes  was  defined  for  each  wall,  floor  or  ceiling. 
An  incident  theta  and  phi  angle.  0,^.^  and  (p^  respectfully  (figure  9).  was  defined  for  each  RAM  lip 
by  the  dot  product  relation  in  equation  (4)  and  the  inverse  tangent  for  the  appropriate  vectors.  A 
scattering  vector  will  then  be  defined  from  the  tip  of  each  RAM  pyramid  to  the  receive  probe.  A 
scattered  theta  and  phi  angle.  0^,  and  0^,  respectfully,  was  defined  in  a  similar  manner  to  the 
incident  angles.  These  angles  will  then  be  inputs  to  the  BISTATIC  subroutine  [13]  in  the  TIPBIST 
file  [13]. 


The  exception  to  the  RAM  patch  method  in  the  high  frequency  model  was  the  back  wall.  The 
incident  angles  are  always  near  the  normal  direction  in  the  chamber  configuration  for  the  cylindrical 
scan.  The  bistatic  scattering  model  will  not  be  appropriate  for  this  wall  so  an  alternative  method 
was  used.  A  specular  incident  and  scattering  vector  was  defined  at  the  conducting  wall  for  each 
turntable  angle  similar  to  the  low  frequency  model.  However,  manufacturer  supplied  data  for  the 
RAM  was  used  for  the  reflection  coefficient. 

B.  RAM  Models  As  stated  above,  the  chamber  verification  was  in  a  chamber  with  24  inch  pyramid 
absorber  on  the  walls,  and  ceiling.  The  exception  was  the  floor  where  48  inch  pyramidal  absorber 
was  placed.  Each  pyramid  absorber  has  a  uniform  base  section  and  the  floor  has  an  intermediate 
layer  of  concrete  that  was  seven  inches  thick  before  the  conducting  shield  was  encountered.  The 
floor  RAM  type  was  a  twisted  pyramid  RAM  which  causes  the  coordinate  system  at  the  tips  to  be 
shifted  by  45  degrees. 

The  geometry  of  the  low  frequency  RAM  model  is  shown  in  figures  10  and  11.  The  total  height  of 
the  RAM  is  L  It  consists  of  a  tapered,  pyramidal  section  of  length  LB  and  a  uniform  section  at  the 
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base  of  thickness  L  •  LB.  The  layer  of  concrete  shown  in  figure  10  was  unique  to  the  floor 
multipath,  otherwise,  the  other  walls,  and  ceiling  will  have  the  pyramid  connected  directly  to  the 
conductor.  The  spacing  between  RAM  tips  was  a,  and  the  width  of  the  pyramid  was  a  function  of 
height,  d(2). 

The  geometry  of  the  high  frequency  model  for  typical  pyramid  RAM  tip  was  shown  in  figure  9.  The 
details  of  how  the  geometry  of  the  pyramid  RAM  was  adapted  to  the  UTD  for  dielectric  corners,  is 
given  in  [13]  and  the  in-depth  discussions  are  left  for  the  interested  reader.  However,  a  brief 
description  will  be  given  as  follows:  A  local  coordinate  system  for  each  edge  was  defined  with  a 
corner  for  each  edge  at  the  apex  of  the  pyramid.  By  the  principle  of  superposition  [28],  the  electric 
field  contribution  from  each  corner  diffraction  can  be  added  to  obtain  the  total  scattering  for  each 
tip.  The  angle,  a,  of  the  RAM  face  is  entered  into  the  subroutine  GEOMETRY  in  the  TIPBIST 
subroutine  which  computes  the  local  coordinate  system  of  each  of  the  four  edges  and  the  unit 
vectors  needed  to  define  the  wedge/corner  diffraction  coordinates,  the  local  polarization  transform 
arxJ  the  dielectric  reflection/transmission  coefficients.  The  unit  vectors  are  transferred  to  BISTATIC 
through  common  variable  statements.  The  exception  to  the  geometry,  as  mentioned  previously, 
was  the  twisted,  pyramid  RAM  as  shown  in  figure  12,  where  the  upper  portion  of  the  pyramid  was 
oriented  at  0,  =  45°  to  the  base  where  the  tip  coordinate  system  was  defined. 

C.  Geometry  of  the  Walls.  Floor  and  Ceiling  The  geometry  of  every  wall,  floor  and  ceiling  follows 
a  similar  format,  so  the  floor  geometry  was  derived  in  detail  and  the  figures  for  the  other  walls,  and 
ceiling  can  be  used  to  visualize  the  relationships  for  that  location.  The  floor  geometry  is  shown  in 
figure  13.  The  source  antenna  image  is  located  by  equation  (7)  and  the  vector  from  the  image  to 
the  probe  is  given  by  equation  (8).  The  specular  angle  at  the  base  of  the  RAM  is  given  in  equation 
(9).  The  vector  to  the  base  of  the  RAM  pyramid  is  found  in  equation  (10)  which  makes  use  of  the 
tangent  formula  in  equations  (11)  and  (12). 
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Figure  9  High  Frequency  RAM  Tip 
Coordinate  System 


Figure  10  RAM  Dimensions  for  Low 
Frequency  Model 


Figure  11  Taper  Geometry  (or  Low 

Frequency  Model  Figure  12  Twisted  RAM  Geometry 


12 


Rjp  -  *  Ay  9  -  A,  S 


(7) 


6^  =  cos' 


* 
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»  +  pL 

(10) 

IF  ^'F 

^F,  _ 

(11) 

Ry  RyJ 
tan  «y  =  -^^  = 


where  a,  is  the  angle  between  the  vector  in  equation  (10)  when  projected  on  the  y  =  0  plane  and 
the  negat'ive  x-axis,  and  is  the  angle  between  the  vector  in  equation  (10)  when  projected  on  the 
X  =  0  plane  and  the  y-axis. 

The  vector  from  the  transmit  horn  to  the  RAM  tip  is  given  by  equation  (13). 

Rrrr  =  Rtf  '  ^  *  yk 


Rtp  =  Rfijr  Rsf  -  Lp  2  R^f 


and  the  vector  from  the  RAM  tip  to  the  receive  probe  is  given  by  equation  (15). 

R^ps  =  f  -  Rtf  *  ^  yk  y 


Now.  the  incident  and  scattered  angles  at  the  tips  are  found  from  equations  (16)  through  (1 9).  which 
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stated  above,  the  figures  for  the  other  walls,  and  ceiling  can  be  used  to  visualize  the  similar 
-eiationships  tar  the  other  walls,  and  ceiling.  See  figure  14  for  the  east  wall,  figure  15  for  the  west 
^3ll,  figure  16  for  the  ceiling  and  figure  17  for  the  back  wall.  The  figures  show  the  appropriate 
coordinate  sysiem  transformation  relative  for  each  wall  and  ceiling.  The  vector  notation  is  consistent 
■nroughout  sc  that  a  change  in  subscripts  is  sufficient  to  derive  the  entire  set. 


“"ne  next  sections  will  describe  the  low  and  high  frequency  simulations  in  detail.  The  flow  of  the 
stectromagneiic  propagation  will  be  described  along  with  more  detailed  description  of  the  RAM 
-modeling  maihnematics.  Computer  resources,  file  management  and  polarization  effects  will  also  be 
discussed. 
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Figure  13  Floor  Geometry 


Figure  14  East  Wall  Geometry 


Figure  15  West  Wall  Geometry 


Figure  16  Ceiling  Geometry 
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Figure  17  Back  Wall  Geometry 
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2.  LOW  FREQUENCY  SIMULATION 

The  low  frequency  model  consists  of  mathematical  descriptions  of  the  electromagnetic  propagation 
in  the  chamber.  An  assumption  was  made  that  the  dominant  scatterer  was  the  specular  reflection 
from  the  walls  so  that  multiple  reflections  and  other  objects  in  the  chamber  will  not  be  included  in 
the  simulation. 

The  flow  of  the  propagating  wave  occurs  in  the  following  manner.  A  stationary,  dipole  antenna 
transmits  to  a  scanning  receive  dipole,  probe  antenna  in  the  anechoic  chamber  and  the  specular 
reflections  from  the  walls,  floor  and  ceiling  of  the  chamber.  The  transmitted  waveform  has  a 
spherical  wave  spreading  and  phase  factor  given  by  the  Green's  function  in  equation  (20), 

(20) 

^  R 


where  R  is  the  distance  from  the  range  antenna  and  k  is  the  free  space  wave  number,  2n/A.  A  ray 
emanates  from  the  range  antenna  directly  to  the  receive  probe  and  also  from  the  range  antenna  in 
the  direction  of  the  specular  reflection  on  the  walls,  floor  and  ceiling  for  each  0  angle.  Each  ray 
was  multiplied  by  the  gain  and  normalized  pattern  factor  of  the  range  antenna.  The  distance  of  the 
specular  paths  was  the  distance  from  the  image  source  to  the  receive  probe.  Each  specular  field 
component  was  multiplied  by  an  attenuation  factor  or  reflection  coefficient  caused  by  the  RAM  on 
the  wall.  The  angle  of  incidence  and  polarization  with  respect  to  the  plane  of  incidence  of  the  ray 
must  be  considered  to  compute  each  reflection  coefficient.  The  E-field  phasors  are  then  summed 
at  the  scan  location  and  converted  to  power  density  by  the  relationship  in  equation  (21), 

5  =  Jll!  (21) 

n 

where  n  is  the  free  space  impedance  of  120r7  ohms,  and  then  plotted  in  dBm/(sq  ft)  versus  scan 
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angle.  The  angle  of  incidence  for  each  ray  to  the  receive  probe  was  then  used  to  compute  the 
receive  antenna  pattern.  Each  ray  E-field  was  muitiplied  by  the  receive  pattern  factor  and  again 
summed,  converted  to  power  density,  and  then  to  received  power  by  the  relationships  in  equation 
(22). 


P,  =  SA 


eff 


4n 


(22) 


where  is  the  effective  aperture  area.  A  is  the  free  space  wavelength,  and  G,  is  the  receive 
antenna  gain.  The  received  power  is  then  plotted  in  dBm  versus  scan  angle. 

A.  RAM  Model  The  low  frequency  RAM  characteristics  viz.  the  bistatic  reflection  coefficient  used 
for  the  specular  ray  attenuation  was  calculated  as  follows:  Since  the  length  of  a  pyramid  was  small 
in  terms  of  wavelength  at  the  low  frequencies,  the  RAM  may  be  modeled  as  a  rough  layer  of 
dielectric  material  with  complex  permittivity  with  periodic  perturbations.  The  layer  of  RAM  will  have 
a  tapered  characteristic  impedance  much  like  a  tapered  impedance  TEM  transmission  line  with  a 
load  impedance,  a  short  circuit  in  most  instances.  The  details  of  the  derivation  of  the  reflection 
coefficient  can  be  found  in  [23]-(25],  and  will  be  described  briefly  here  without  proof.  The 
constituent  properties  of  the  material  are  modeled  as  an  inhomogeneous,  uniaxially  anisotropic 
absorbing  dielectric.  The  dielectric  and  permittivity  constants  will  be  in  a  tensor  format  as  in  matrix 
equations  (23)  and  (24). 


(24) 


P  = 


0 

0 


0 

0 


0 

0 


The  RAM  pyramids  may  be  modeled  as  a  periodic  array  of  rods  of  spacing  a  (figure  18)  that 
become  tapered  to  form  a  pyramid.  The  relative  volume  of  space  occupied  by  the  pyramid  at  a 
cross  section  in  the  z  plane  is  given  by  g  =  d(z)/a.  When  a  becomes  small  in  terms  of 
wavelengths,  the  material  properties  can  be  modeled  by  equations  (25)  through  (28), 


e.  = 
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ey  = 


(i-g^) 
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l+gr2. 


+  g% 

\^o  * 

2  (ga-ep) 

(l+g^)e^  +  (l+gr-) 


(25) 

(26) 


(27) 


1+g^ 
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(i*g^)  \io  *  (1-9^^) 


(28) 


where  e,  and  are  the  relative  permittivity  and  permeability  parameters,  respectively,  and  and 
are  the  permittivity  and  permeability  of  free  space,  respectively.  It  is  assumed  that  all  of  the 
materials  used  are  non-ferrous,  so  //,  =  For  an  E-field  that  is  perpendicular  to  the  plane  of 
incidence  the  effective  material  properties  are  given  in  equations  (29)  and  (30), 


'eff 


ey  - 


(29) 


(30) 


and  in  equations  (31)  and  (32)  for  an  E-field  that  is  parallel  to  the  plane  of  incidence. 
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(31) 


e^^i^sin^e 


(32) 


The  equivalent  TEM  transmission  line  parameters  of  the  characteristic  impedance,  IJ,z),  and  the 
propagation  constant,  y{z),  are  given  by  equations  (33)  and  (34). 


N 


(33) 


y  {z)  -  J  (2)  e.gff  (  z) 


It  was  assumed  that  the  reflection  coefficient  was  small  enough  that  smaller  order  terms  can  be 
neglected  [32],  the  reflection  coefficient  from  a  typical  wall  is  given  by  equation  (35), 

r(o)  =  ri(o)  +  r2(o)  (^s) 


where 


ri(o)  = 


(36) 


r2(0)  =  f^'^N(z')  dz' 

J  0 


(37) 


Niz)  = 


2Z^(z) 


(38) 


and  rL  =  -7.  Since  the  material  was  continuous  between  the  pyramid  tapered  layer  and  the 
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uniform  bottom  layer,  the  reflection  coefficient  was  assumed  to  be  zero. 


The  floor  has  an  intermediate  concrete  layer  with  thickness  d^on  and  was  modeled  by  equation  (39), 


r^(o)  =r^(o)  +rp(o)  +rp(o)  (39) 


where. 


r^KO) 


(40) 
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con 


^  con  ^  Ahs 
^  con  ^  ^  ahs 


(41) 
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rp(o)  = 


y  2  ~  J  ^'J^con^con 


(42) 

(43) 

(44) 

(45) 


multiple  reflections  are  neglected  between  layers  since  the  reflections  are  small,  and  ticon^Hc 
€^^=10€g  are  assumed  as  the  constituent  properties  of  the  concrete  layer  (33],  where  the 
permittivity  has  been  estimated  based  on  the  permittivities  of  the  concrete  materials.  The  reflection 
between  the  pyramidal  and  uniform  layer  are  again  assumed  to  be  zero. 


B.  Computer  Resources  The  simulation  has  been  developed  on  a  Z486  PC  and  uses  a  MATHCAD^ 
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version  3.0  software  package  [34],  Since  a  single  file  will  contain  a  maximum  of  50  variables,  which 
decreases  as  the  number  of  iterations  or  integrals  increases,  several  files  were  cascaded  to 
complete  the  computation.  The  files  were  connected  by  WRITE  and  READ  statements  which  passed 
the  relevant  arrays  between  files.  A  total  of  26  files  were  used  in  the  low  frequency  model,  and  the 
flow  chart  of  figure  19  shows  the  order  of  the  files.  A  check  was  added  to  the  beginning  of  each 
file  to  assure  the  order  was  maintained  and  a  statement  of  where  to  go  next  was  placed  at  the  end 
of  each  file  for  convenience. 

The  integrals  evaluated  in  the  RAM  model  in  MATHCAD®  use  a  numerical  algorithm  called  a 
Romberg  integration.  The  lower  limit  of  the  reflection  coefficient  integral  had  to  be  offset  slightly 
from  zero  to  avoid  a  convergence  problem  with  the  algorithm. 

Q  Polarization  and  Orientation  Considerations  The  polarization  for  the  low  frequency  simulation 
was  vertical  to  the  chamber  floor.  The  propagating  E-field  was  represented  in  polar  coordinates  as 
Eq.  The  floor  polarization  remains  the  same  throughout.  The  ceiling  orientation  of  the  E-field  was  - 
Ee  when  referenced  to  the  local  coordinate  system.  The  east  and  west  wall  E-field  orientation  were  - 
E,  and  E,  respectfully  for  the  local  coordinate  system.  The  east,  west  and  back  walls  have  a  180 
degree,  boundary  value,  phase  shift  in  the  low  frequency  model  consistent  with  image  theory. 

The  reflection  coefficient  for  each  polarization  component  was  referenced  to  the  plane  of  incidence. 
The  floor  and  ceiling  E-field  consistently  has  only  a  component  that  was  parallel  to  the  plane  of 
incidence  (figure  20).  However,  the  east,  west  and  back  wall  each  have  a  parallel  and  perpendicular 
component.  The  derivation  for  the  components  of  the  east  wall  will  be  shown  here.  The 
components  for  the  other  walls  were  derived  and  used  in  a  similar  manner.  Figure  21  shows  the 
field  orientation  and  the  plane  of  incidence  referenced  to  the  east  wall.  The  unit  vectors  along  the 
direction  of  propagation  with  respect  to  the  plane  of  incidence  are  given  by  equations  (46)  and  (47). 
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Figure  18  Transverse,  Periodic,  Rectangular, 
Lossy  Rods  Used  in  Low  Frequency 
Formulation 


Figure  19  Flow  Chart  for  Low  Frequency  Model 
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Figure  20  Floor  Polarization  Orientation 


Figure  21  East  Wall  Polarization  Orientation 
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The  E-field  was  converted  to 


rectangular  coordinates  by  multiplication  by  the  vector  in  equation  (48). 
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The  parallel  and 


perpendicuiar  components  are  then  given  by  equations  (49)  and  (50). 
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(49) 


E,  =  E^'U^^ 
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The  parallel  and  perpendicular  fields  are 
then  summed. 


then  multiplied  by  the  appropriate  rcBcction  coellicient  and 


D.  n. . ..Ut  nene-ive  Antennas  A  half-wave  dipole  oriented  vertically  or  on  the  z-axis  ol  the 

turntable  coordinate  system,  was  used  forthe  source  and  recewe  antennas.  A  hall-wave  dipole  has 
a  directivity  of  t.64  or  2.1  dBi  (27,28]  and  has  an  omnielirectional  pattern  in  p.  Each  ray  in  the 
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simulation  was  multiplied  by  a  transmit  normalized  antenna  pattern  defined  by  the  angle  0  of  the 
ray.  The  receive  ray  was  multiplied  by  the  normalized  dipole  pattern  if  the  received  power  was  the 
displayed  quantity.  The  elevation  pattern  for  a  typical  half-wave  dipole  is  given  by  equation  (51). 

cos(-|cosej  (51) 

dipole  sin0 

The  derivation  of  the  direct  path  and  floor  angles  0,  and  0,  is  shown  here.  Equations  for  these 
angle  for  the  other  walls  were  derived  in  a  similar  manner.  The  angles  for  the  direct  path  are  given 
by  equations  (52)  and  (53), 


=  cos-1 
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and  the  angles  for  the  floor  are  given  by  equations  (54)  and  (55). 
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E.  Ntimher  of  Data  Points  The  number  of  data  points  in  a  scan  was  determined  by  the  Nyquist 
criteria  (26)  [35],  and  also  by  the  FFT  algorithm  in  the  cylindrical  wave  expansion.  For  a  cylindrical 
scan  the  maximum  angle  increment  determined  by  the  Nyquist  criteria  is  given  in  equation  (56), 
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which  for  a  2ij  radian  scan  reduces  to  kR  ^  N^2.  The  FFT  algorithm  requires  that  the  number 
of  data  points  by  a  factor  oi  N  =  2^,  where  n  is  an  integer. 

The  high  frequency  model  will  be  described  next.  The  frequency  range  of  the  high  frequency  model 
will  be  applied  was  500  MHz  to  18  GHz.  This  model  can  also  be  applied  to  the  low  frequencies 
since  the  model  assumes  a  pyramid  of  infinite  length;  but,  since  the  specular  reflection  was  so  large 
at  the  lower  frequencies,  the  diffraction  effects  became  less  significant. 


3.  HIGH  FREQUENCY  SIMULATION 

The  high  frequency  model  was  also  a  mathematical  description  of  the  propagation  and 
scattering/multipath  interference  in  a  rectangular  anechoic  chamber.  The  overall  format  was  the 
same  as  the  low  frequency  model.  The  major  differences  are  the  way  the  RAM  was  modeled  and 
the  type  of  range  antenna  used.  The  RAM  model  assumes  that  the  RAM  tips  are  the  dominant 
sources  of  scattering.  Since  there  are  many  RAM  tips,  the  number  of  GO  rays  that  are  modeled  and 
summed  are  significantly  increased  compared  to  the  low  frequency  model.  For  instance,  a  20  x  20 
RAM  patch  has  1602  rays  versus  a  maximum  of  six  rays  in  the  low  frequency  model. 

The  power  flow  in  the  high  frequency  model  was  similar  to  the  low  frequency  model  except  that 
image  theory  was  not  applied.  This  means  that  the  phase  shift  will  not  occur  since  the  scattering 
was  coming  from  the  tips  and  not  the  conducting  wall.  Also,  the  distances  used  in  the  spherical 
wave  spreading  attenuation  was  calculated  from  the  vector  from  the  range  antenna  to  the  RAM  tips 
plus  the  vector  from  the  RAM  tips  to  the  receive  antenna.  An  assumption  was  also  made  that  these 
two  distance  vectors  are  nearly  equal  when  computing  the  reflection  coefficient.  The  spherical  wave 
spreading  was  added  in  this  model  because  it  was  not  included  in  the  TIPBIST  formulation  in  [13]. 
The  diffraction  model  assumes  that  the  incident  and  reflections  are  far  enough  away  that  a  plane 
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wave  assumption  was  made.  Since  only  the  relative  value  between  the  incident  and  scattered 
energy  was  used  to  compute  the  reflection  coefficient,  the  input  E-field  was  set  to  unity  with  the 
appropriate  polarization  and  the  reflection  coefficient  was  computed  from  E*^/E'"‘. 


A.  RAM  Model  The  input  to  the  TIPBIST  GEOMETRY  subroutine  wedge  angle  and  the  inputs  to 
the  BISTATIC  subroutine  are  Ee~  e;"  0'"*.  0^  e*^.  and  0*“'.  These  quantities  were  computed 
in  the  Geometry  section  of  this  report.  The  RAM  model  was  a  summation  of  four  dielectric  edges 
with  a  corner  diffraction  term  for  each  edge  at  the  apex  of  the  pyramid.  The  four  edges  are 
modeled  separately  with  a  three  dimensional,  dielectric,  UTD  algorithm  with  the  comer  diffraction 
term  added.  The  details  of  this  algorithm  are  developed  shown  in  reference  [13],  where  the 
development  of  the  various  stages  of  the  model  was  found.  An  overview  of  the  relevant  algorithm 
used  in  TIPBIST  will  be  shown  below  with  the  same  notation  as  that  of  reference  [13].  The  overall 
equation  will  be  shown  for  one  edge  of  the  RAM  and  each  component  of  the  equation  will  be 
explained.  Since  a  great  deal  of  the  formulation  was  spent  in  transforms  between  coordinate 
systems  and  in  model  interfaces  (from  models  that  were  developed  separately),  the  derivation 
shown  here  will  be  brief. 

The  overall  equation  for  the  scattered  E-field  in  BISTATIC  is  given  in  equation  (57). 
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Poc 
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where, 
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F(x)  =  2jv/x  e'^’’dt 

J  y/x 


(66) 


a*  (p) 


=  2cos^ 


-  p 

2 


atxJ  N*  are  integers  that  are  derived  by  rounding  off  the  solutions  of  equation  (68), 


(67) 


InnN*  -  P  =  ±7t 

n  =  2-WA/n  (^5) 

Also,  WA  is  the  wedge  angle,  k  is  the  wavenumber,  e,  is  the  free  space  permittivity,  is  the 
complex  permittivity  of  the  RAM  material.  The  o  subscript  refers  to  the  incident  face  of  a  wedge  and 
n  refers  to  the  shadow  face  of  the  wedge.  The  function  F(x)  is  the  transition  function  used  in  the 
standard  UTD  formulation;  however,  the  ratio  L/L,  is  considered  constant  in  BISTATIC.  The  C. 
through  are  the  corner  diffraction  terms  for  four  wedges,  and,  the  R|  and  are  the  reflection 
coefficient  terms  used  to  account  for  the  dielectric  boundaries  of  the  wedge.  The  parameters  S,  Bo,., 
Sj,  and  So  are  defined  in  figures  23  through  25. 


B.  Back  Wall  As  mentioned  in  a  previous  section,  the  incident  angle  to  the  back  wall  was  assumed 
to  be  always  near  normal  incidence  and  was  therefore  modeled  by  manufacturer  supplied  data.  The 
back  wall  attenuation  for  the  RAM  was  specified  as  -50  dB.  The  free  space  attenuation  for  the  back 
wall  was  significantly  larger  than  the  other  walls,  so  the  assumption  will  not  affect  the  outcome  of 
the  results. 


C.  Computer  Resources  The  high  frequency  model  was  programmed  in  Fortran  on  a  VAX 
mainframe  computer  and  linked  with  the  TIPBIST  file.  A  flow  chart  of  the  program  is  shown  in  figure 
26.  The  amount  of  memory  needed  for  the  large  arrays  from  the  RAM  patches  necessitated  a 
mainframe  computer. 
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Figure  26  Flow  Chart  for  High  Frequency  Model 
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D.  Receive  and  Source  Antenna  The  receive  probe  antenna  was  a  half-wave  dipole  as  in  the  low 
frequency  model.  The  source  antenna  for  the  high  frequency  model  was  a  TE„  diagonal  horn.  The 
mathematical  model  for  the  radiation  pattern  was  developed  by  Love  [37].  The  radiation  patten  was 
symmetric  in  the  principal  planes  though  somewhat  degraded  in  the  intercardinal  planes.  The 
directivity  of  the  horn  was  found  with  the  equation  D=4nA,/A  with  an  efficiency  of  0.81 .  The  mode 
structure  in  the  aperture  is  shown  in  figure  27  and  the  local  coordinate  system  is  shown  in  figure 
28.  A  coordinate  transformation  was  necessary  for  each  wall  to  the  local  coordinate  system  of  the 
wall.  The  transformation  will  be  shown  for  the  direct  path  and  the  floor,  whereas  the  formulation  for 
the  other  walls  was  similar. 


The  E-field  for  the  TEq,  and  TE,o  modes  are  orthogonal/independent  so  that  they  may  be  calculated 
separately.  The  resultant  field  was  then  found  from  the  combination.  The  E-field  for  each  mode  is 
given  by  equations  (70)  and  (71). 
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where, 

u  =  sin  e  (72) 

A 
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and  d  is  the  length  of  the  side  of  the  square  aperture.  The  coordinate  transfer  for  the  direct  path 
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is  given  in  equations  (74)  through  (76). 
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The  coordinate  transfer  for  the  floor  path  is  given  in  equations  (77)  through  (79). 


^fDH  ~ 
^fDH  ~ 


cos'^ 

cos"^ 


i? 


fD/f, 


t  l-^rfrl ) 


'■fDH, 


\  l-^rrrl  sin  Q/xjh/ 


^  ®  fDH ' 

Ef  =  0 
■f* 

^fDH,  ~  ~^TFTji 

^fDHy  ~  ^^TFTy  *  ^TFt)  ^ 

^fDM^  ~  ^^TFTf  ~  ^Trr)  ! 


(77) 


(78) 


(79) 


36 


E.  Number  of  Data  Points  The  number  of  data  points  used  the  same  criteria  as  the  low  frequency 
model. 


4.  MODAL  EXPANSION  PROCESSING 

A  modal  expansion  of  the  E-field  can  be  used  to  determine  the  fields  at  any  other  radius  than  the 
measurement  radius.  The  fields  converge  in  phase  when  a  radiating  source  was  encountered,  and 
can  be  used  to  locate  scattering  sources  in  the  chamber.  The  methodology  of  the  simulation  and 
verification  was  to  compare  the  predicted  data  to  measured  data.  When  a  discrepancy  was 
experienced,  a  modal  expansion  was  performed  on  the  measured  data  to  determine  where  the 
scattering  was  coming  from.  It  was  hopeful  that  scattering  sources  other  than  the  RAM  covered 
wails  can  be  located  and  corrected. 


The  modal  expansion  used  was  found  in  reference  [26]  and  is  given  by  equation  (80), 


FFT-'^ 


FFT[E,(R,<i>)]  g  ^,2) 

t  {kR) 


(80) 


where  R  is  the  scan  radius,  E2(R,p)  is  the  measured  data,  FFT  is  a  Fast  Fourier  Transform,  and 
is  a  Hankel  function  of  the  second  kind.  The  distance  resolution  of  the  measurement  was  obtained 
by  sweeping  several  frequencies  around  the  center  frequency.  After  the  modal  expansion,  the 
frequency  components  of  the  electric  field  are  added.  The  phasors  of  the  field  add  when  a  source 
was  encountered  in  the  chamber  and  create  a  maximum  in  the  field  magnitude.  When  there  were 
no  prevalent  scattering  sources,  the  electric  fields  add  in  phase  randomly  and  do  not  typically  create 
any  significant  magnitude. 
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5.  MEASUREMENT  CONFIGURATION 


The  measurements  were  conducted  in  a  rectangular  chamber  as  described  in  the  geometry  section. 
An  automatic  vector  network  analyzer  (Hewlett  Packard  HP8510)  coherent  detection  system  was 
used  to  receive  and  collect  the  data.  A  block  diagram  of  the  equipment  is  shown  in  figure  29.  A 
Flam  and  Russell  959  system  was  used  to  process  the  data  for  an  antenna  pattern  format  and 
controlled  the  antenna  positioning  system.  The  80  foot  chamber  turntable  was  not  available  during 
the  measurements,  so  a  Scientific  Atlanta  (SA)  azimuth  positioner  was  substituted  for  the 
measurement.  The  location  of  the  positioner  in  the  chamber  is  shown  in  figure  30.  An  arm 
constructed  of  low  density  wood  was  fabricated  to  provide  the  cylindrical  scan  (figure  31).  The 
positioning  equipment,  wooden  arm  and  input  cable  was  covered  with  available  RAM  so  that 
additional  scattering  would  not  interfere  with  the  measurements.  The  RAM  was  placed  to  cover  the 
input  cable  to  the  antennas  as  much  as  possible  to  prevent  scattering  from  the  cable  from 
interfering  with  the  measurements. 

The  antennas  used  for  the  simulation  verification  measurements  were  dipoles.  A  sleeve  dipole  tuned 
to  10  GHz  was  used  for  the  high  frequency  measurements.  The  sleeve  dipole  was  constructed  of 
0.141  coaxial  cable  with  a  brass  sleeve  choke  to  form  the  lower  half  of  the  dipole  and  suppress 
input  cable  currents.  The  dipoles  for  the  low  frequency  measurements  were  from  a  calibrated 
EMCO  Model  31 21 C  Tuned  Dipole  Antenna  set.  The  dipoles  were  tuned  and  measured  (figure  32 
and  33)  with  an  HP8757  Scalar  Network  Analyzer. 

The  antenna  used  for  the  diagnostic  measurements  was  a  broadbeam  horn.  A  TECOM  Type  No. 
201 187-4  Dual-Polarized,  Quad-Ridge  Horn  was  selected.  The  cross  polarized  port  of  the  horn  was 
terminated  with  a  50  ohm  load.  The  beamwidth  of  the  horn  was  80°  to  40°,  which  was  similar  to 
an  open-ended  waveguide.  The  horn  antenna  was  substituted  as  close  as  possible  to  the  exact 
location  of  the  dipole  at  the  top  of  the  positioning  arm.  The  horn  fixture  below  the  dipole  location 
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Figure  29  Measurement  Equipment  Block  Diagram 
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Figure  30  Location  of  the  S/A  Positioner  in  the  Chamber 


Figure  31  Positioner  Configuration  with  Antenna  Fixture 
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shown  in  figure  30  W3S  fabricsted  to  hold  heavier  and  more  bulky  horn  antennas  for  facility 
measurements  not  related  to  this  project. 

The  HP8510  is  a  coherent,  ratio  measurement,  though  absolute  power  measurements  are  a 
capability  of  the  system.  The  simulated  plots  were  corrected  for  power  levels;  however,  the  actual 
measurements  were  adjusted  in  scale  to  the  power  level  of  the  simulation  so  that  the  measurements 
can  be  overlaid  for  comparison. 

The  next  section  is  a  comparison  of  the  predicted  and  measured  results.  An  exception  will  be  noted 
concerning  two  antennas  in  the  chamber  that  were  present  during  the  measurement. 

6.  DATA  Comparisons 

An  exception  to  the  simulation  model  was  made  to  account  for  two  large  antennas  that  were 
present  in  the  chamber  that  could  not  be  removed  for  the  measurements.  The  antennas  were  Flam 
and  Russell  6418,  diagonal  horn  antennas  with  apertures  of  2.5  ft.  side  dimensions.  The  location 
of  the  first  horn  was  at  a  radius  from  the  center  of  the  turntable  of  134.5  ft.  at  0  =  123.88°.  The 
location  of  the  second  horn  was  at  a  radius  of  138.7  ft.  and  0  =  -124.34°.  An  approximation  of  the 
scattering  of  the  antennas  found  in  references  [38]  and  [39]  was  added  to  the  simulation  model  with 
an  antenna  cross  section  approximation  given  in  equation  (81). 

o  =  Gain  •  Area  = 

A 


A  bistatic  radar  range  equation  was  used  to  model  the  antennas  where  the  distance  from  the 
antennas  to  the  receive  probe  antenna  changed  as  the  scan  occurred. 

An  overlay  of  the  low  frequency,  300  MHz,  simulated  and  measured  data  is  shown  in  figure  34.  The 
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direct  path  was  also  plotted  to  show  the  effect  of  free-space  attenuation  on  the  measurement.  The 
standing  wave  for  this  chamber  configuration  was  approximately  6  to  10  dB. 

A  second  low  frequency  simulation  was  conducted  at  300  MHz  without  the  chamber  antennas  as 
shown  in  figure  35.  The  floor  reflection  was  the  most  significant  scatterer  of  2.8  dB.  The  ceiling  was 
expected  to  contribute  1.1  dB.  The  other  walls  contribute  much  less  in  magnitude  but  at  a  higher 
spatial  frequency  since  the  path  length  difference  between  the  direct  and  secondary  path  was  much 
greater  for  the  other  walls.  The  path  length  difference  for  the  floor  was  approximately  one 
wavelength.  The  change  in  the  path  length  difference  for  the  entire  floor  scan  was  only 
approximately  3  inches.  Because  of  these  effects  and  the  floor  being  the  most  significant 
contributor,  the  simulated  data  always  stays  positive  for  the  entire  scan. 

The  high  frequency,  10  GHz.  prediction  and  scan  are  shown  overlaid  in  figure  36,  and  separately 
in  figures  37  and  38,  respectfully.  The  simulated  data  was  limited  to  a  21  x  21  patch  of  RAM  tips 
due  to  the  limitations  of  the  VAX  disk  pack.  The  ripple  levels  were  between  1  and  2  dB  peak-to- 
peak  in  the  predicted  data. 

The  modal  expansion  data  was  expanded  to  several  radii  at  the  locations  where  scattering  sources 
in  the  chamber  are  known  to  be.  The  number  of  data  points  used  was  512.  These  included  the 
corners  of  the  chamber,  the  chamber  antennas  at  the  opposite  ends  of  the  chamber,  and  the 
turntable  edges.  Figure  39  shows  the  only  plot  that  revealed  any  significant  information  in  the  -90° 
to  -180°  sector.  This  scattering  source  was  probably  the  chamber  antennas  mentioned  previously 
and  included  in  the  low  frequency  simulation, 

A.  Data  Correlation  The  levels  of  the  ripple  agrees  well  between  the  simulation  and  the 
measurements  when  the  chamber  antennas  were  included  in  the  simulation  for  approximately  90 
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Figure  36  High  Frequency  Comparison  of  Simulated  and  Measured  Data 
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%  of  the  scan.  However,  measurements  without  the  antennas  were  not  possible  in  the  chamber  to 
verify  the  simulation  without  these  antennas. 


The  high  frequency  measurements  show  approximately  a  1  dB  higher  peak-to-peak  ripple  than  the 
simulation  in  80%  of  the  scan.  In  the  ±  40°  to  70°  sector  a  slightly  larger  standing  wave  of 
approximately  4  dB. 

B.  Rpsolvina  Di.screpancies  The  discrepancies  in  the  high  frequency  measurements  were  possibly 
caused  by  the  manufacturing  process  of  the  twisted  pyramidal  RAM  on  the  floor.  It  was  determined 
that  it  was  difficult  to  use  the  modal  expansion  diagnostic  method  in  the  direction  of  the  source 
antenna.  It  was  not  possible  during  the  measurements  to  replace  the  twisted  pyramid  RAM  with 
other  types  of  RAM  to  determine  the  source  of  the  low  frequency  ripple.  The  modal  expansion  was 
attempted  to  resolve  discrepancies  in  the  high  frequency  model,  although  none  were  found.  The 
resolution  of  the  modal  expansion  was  possibly  affected  by  the  1  GHz  frequency  sampling  interval. 


7.  CONCLUSIONS 

This  project  was  an  attempt  to  improve  the  quality  of  chamber  characterizations  in  an  anechoic 
chamber.  Computer  simulated  data  was  compared  to  measured  data.  Good  agreement  was  found 
between  the  simulation  and  measurement  for  the  low  frequency  data  when  the  Flam  and  Russell 
6418  antennas  were  included  in  the  simulation.  There  was  good  agreement  in  the  high  frequency 
data,  except  for  a  slight  variation  from  a  low  frequency  standing  wave. 

A  Recommendations  The  following  recommendations  are  made  that  need  to  be  investigated  to 
improve  the  simulation  and  measured  data: 
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1.  The  physical  characteristics  of  the  chamber  such  as  the  chamber  corners,  chamber 
antennas,  input  cables,  the  chamber  turntable  edges,  and  the  wall  behind  the  range 
antennas  can  be  included  in  the  modeling. 

2.  Better  constituent  parameters  need  to  be  obtained  for  the  RAM  and  concrete  materials  used 
in  the  simulations. 

3.  A  much  finer  sampling  interval  was  needed  for  the  modal  expansion  resolution. 

4.  The  low  frequency  files  need  to  be  combined  into  one  program. 

5.  Linking  of  more  than  one  disk  pack  on  the  VAX  computer  would  allow  a  larger  patch  of 
RAM  in  the  high  frequency  simulation. 

6.  The  modal  expansion  scan  data  would  be  much  more  informative  in  a  waterfall  plot  format. 

7.  The  Flam  and  Russell  641 8  antennas  need  to  be  removed  from  the  chamber  when  they  are 
not  being  used. 

B.  Summan/  A  simulation  of  an  anechoic  chamber  with  a  cylindrical  scan  was  developed.  Ray 
tracing  was  used  to  model  the  wall  scattering  with  a  RAM  model  used  to  predict  the  wall  scattering. 
The  simulation  was  broken  up  into  a  low  and  high  frequency  model  due  to  the  frequency  dependant 
nature  of  the  RAM  scattering.  The  low  frequency  simulation  was  developed  using  a  single,  specular 
ray  for  each  wall  and  the  RAM  was  modeled  using  a  layered,  tapered,  transmission  line  method  with 
an  anisotropic,  inhomogeneous  tensor  RAM  model.  The  high  frequency  model  was  developed  using 
multiple  rays  on  each  wall  to  the  tips  of  the  RAM.  The  tips  of  the  RAM  were  modeled  using 
TIPBIST.FOR,  a  bistatic  scattering  code  for  a  single  tip  using  a  dielectric,  corner  UTD  formulation. 
Each  tip  is  a  combination  of  four  edges.  Each  model  sums  all  of  the  vectors  at  the  receive  probe 
antenna.  A  simulation  model  was  included  for  the  range  and  receive  antennas. 

Measurements  were  made  in  an  anechoic  chamber  and  a  validation  of  the  simulated  data  was 
conducted.  The  measured  data  correlated  well  with  the  simulation. 


48 


A  cylindrical  scan,  modal  expansion  was  developed  to  resolve  discrepancies  between  the  simulated 
and  measured  data.  The  measurements  were  conducted  for  incoming  waves  only  using  more 
directional  antennas  than  the  characterization  probes. 
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